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TECHNICAL NOTE 2982 

SUPERSONIC FLOW PAST OSCILLATING AIRFOILS 
INCLUDING NONLINEAR THICKNESS EFFECTS 
By Milton D. Van Dyke 


SUMMARY 


A solution to second order in thickness is derived for harmonically 
oscillating two-dimensional airfoils in supersonic flow. For slow oscil- 
lations of an arbitrary profile, the result is found as a series includ- 
ing the third power of frequency. For arbitrary frequencies, the method 
of solution for any specific profile is indicated, and the explicit 
solution derived for a single wedge. 

Nonlinear thickness effects are found generally to reduce the 
torsional damping, and so enlarge the range of Mach numbers within which 
torsional instability is possible. This destabilizing effect varies only 
slightly with frequency in the range Involved in dynamic stability analy- 
sis, but may reverse to a stabilizing effect at high flutter frequencies. 
Comparison with a previous solution exact in thickness suggests that 
nonlinear effects of higher than second order are practically negligible. 

The analysis utilizes a smoothing technique which replaces the 
actual problem by one involving no kinked streamlines. This strategem 
eliminates all consideration of shock waves from the analysis, yet 
yields the correct solution for problems which actually contain shock 
waves . 


INTRODUCTION 


As linearized supersonic -flow theory is increasingly applied to 
problems of unsteady motion of lifting wings, the results are sometimes 
advanced with the warning that they may be significantly affected by 
nonlinear effects of thickness. Such caution is justified because it is 
known that even for steady flow, linearized theory is often inadequate 
for predicting the pitching moment - and prediction of moments is one of 
the main objectives of unsteady-flow theory. It may be anticipated that 
nonlinear effects will become increasingly important as the Mach number 
falls toward unity, particularly for slow oscillations. 
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In "the present work the effects of thickness are determined for a 
harmonically oscillating two-dimensional airfoil by calculating the 
second-order solution. This is the counterpart for uns teady motion of 
the well-known steady-flow result of Bus emann (ref. l) . First, for 
slow oscillations a solution is found for an airfoil of arbitrary pro- 
file. The result is given as a series which includes terms up to the 
third power of the frequency. Second, for arbitrarily high frequencies 
it is shown that a solution can be found for any specific airfoil, and 
the solution is carried out explicitly for a single wedge. Finally, 
comparison is made with a previous solution for the wedge which is exact 
with respect to thickness (refs. 2 and 3 ), in order to assess the effects 
of nonlinear terms of higher than second order. 

Extensive use is made of a smoothing technique, which replaces the 
actual problem by one having no kinked streamlines. This strategem, 
which has been used previously and may prove useful in future problems, 
eliminates all consideration of shock waves from the analysis. It, 
nevertheless, leads to the correct second-order solution for the actual 
problem, which does involve shock waves. 


METHOD OF ANALYSIS 
Statement of Problem 


Consider a uniform supersonic stream flowing past a Bharp-nosed 
airfoil which executes prescribed harmonic oscillations. We shall be 
concerned with calculating the instantaneous pressure at the surface of 
the airfoil and, hence, the unsteady lift and pitching moment. It is 
usually the custom to consider the oscillation to be compounded of a 
rotation (say about the midchord) and a vertical translation (plunging) , 
which are not generally in phase . However, we shall here limit atten- 
tion to the special case of the airfoil pivoted about an arbitrary point 
downstream of its leading edge. This is no real restriction because the 
plunging case can be recovered by letting the distance to the pivot tend 
to infinity and the angular amplitude of oscillation tend to zero, their 
product remaining finite. 

Although the iteration procedure to be employed yields a formal 
result for any Mach number greater than unity, the solution probably 
breaks down when the flow becomes sonic at any point. Since this occurs 
at a Mach number somewhat higher than that for bow wave detachment, the 
upper and lower surfaces of the airfoil operate independently in the 
probable range of validity of the solution. It is therefore sufficient 
to consider only the half field of flow lying above the airfoil, and 
this viewpoint will be adopted henceforth. 
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Choose the origin of coordi- 
nates at the mean position of the 
leading edge, with the x axis 
extending in the direction of the 
free stream (sketch (a)). Then 
it is convenient to describe the 
upper surface of the airfoil in 
its mean (zero angle of attack) 
position by 



Sketch (a).- Coordinate system for 
oscillating airfoil. 


y = Y(x) = e g(x) (l) 

Here e is a small parameter representative of the airfoil thickness, 
so that the function g is of order unity. Now let the airfoil pivot 
about a point lying a distance b downstream from the leading edge, and 
perform harmonic oscillations of frequency w and amplitude ccq, so 
that the angle of attack, which is the angle between the instantaneous 
position of the airfoil and its original mean position (sketch (a)), is 
given by 


a = 0,3 cos wt = a 0 e lwt (2) 

(Here, as in all that follows, it is implied that actual physical 
quantities are given by the real parts of their complex representations.) 
Then, at any instant the moving upper surface of the airfoil is described 
implicitly, to second order in thickness and angle of attack, by 

y = € g(x-ay) - (x-b)a. (3a) 

Now, if the airfoil is smooth (so that the function g is differenti- 
able), the equation of the surface can be written explicitly as 

y = e g(x) - a^e^ (x-b) (3b) 

with an error of order (e 2 a,scc 2 ), which is of third order and, conse- 
quently, negligible in the present second-order analysis. 


Perturbation Equation 


The entropy changes due to shock waves are of third order in the 
airfoil thickness and angle of attack. Hence, to second order the flow 
is irrotational and isentropic. Because it is irrotational, there exists 
a potential function ft whose gradient yields the velocity vector: 
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?= grad fl, v = ny | M 

Bernoulli's equation for plane unsteady flow can be written (from 
eqs. (l^.O^) and (9-06) of ref. 4) 

fl t + i (u 2 + v 2 ) + — = + I u 2 (5) 

2 7-1 7-1 2 

Here a is the speed of sound, and a Q its value in the free stream, 
where the flow velocity is U. Differentiating this expression with 
respect to time t, and using the fact that d[a 2 /(7-l)] = a 2 dp/p 
(ref. 1*-, eqs. (9*03) and (9-06) gives 


fltt + uut + wt + ^- Pt = 0 (6) 

This, together with the corresponding results obtained by differentiating 
with respect to x and y, can be used to eliminate derivatives of the 
density from the continuity equation (ref. 4, eq. (7-08*2)) 

p t + ( pu) x + ( pv )y = 0 (7) 

The result is that the velocity potential satisfies the equation 

(a 2 -ftx 2 ) fixx + (a 2 -fly 2 )flyy - 2fix;flyflxy “ ^flxflxt - 2nyDyt - fltt = 0 (8a) 

where, from equations (4) and (5) , 

a 2 = ao 2 + (U 2 - flx 2 - 

Now introduce a perturbation potential $>, 
by the free-stream velocity U, by setting 

0 (x,y,t) = U [x + $ (x,y, t) ] (9a) 


fly 2 - 2fl t ) (8b) 

normalized through division 


so that the velocity components are given by 


u 

U 


1 + 0 X 



(9b) 
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Then substituting into equations (8) gives 




(r-i) (« x + hi + (t^+jyy) 


( 2 V $ x 2 ) $ xx + V^YY + 2 ( 1 + $ x) Vxy + 


u ( ®x 4 xt + Vyt } 


( 10 ) 


For purposes of a second-order solution (and to higher order a potential 
does not exist), the triple products on the right-hand side can be dis- 
regarded. Thus, the perturbation equation becomes finally 


V 132 ®**- 2 - $®tt - * 


~) ( ®xs + ^yy ) + 


where p 2 = M 2 - 1. 


[( 7 - 1 )( 4 > x + ~) 

+ 2®y ®xy + ^ (®x®xt + ®y®yt)J(ll) 


Pressure Relation 


Dividing the Bernoulli equation (eq. (5)) by a 2 /( 7 ~l) gives 

2Z± ( x , (12 ) 


a 

a o 




The flow is isentropic to second order, so that 


£-(&) 


7 - 1 , 


Po 


|p 0 U 2 




and it follows that the pressure coefficient at any point is given by 

7 


c_= p - p ° 


s I peY 2 


2 

7 M 2 


1 + 


7-1 


M* (l- 7 " 1 


( 13 ) 
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Substituting the -velocity components of equation (9b) , expanding in 
series, and retaining only squares and products of perturbation quanti- 
ties gives finally 


Cp ■ - .2 1 + -® y 2 + 2lP* x £ + M 2 ( £=) (14) 

Here the second-order solution is required only for evaluating the first 
two terms; the others are given correct to second order by linear theory. 


Smoothing of Problem 


From the leading edge and from any subsequent corner of the airfoil 
spring shock waves or Prandtl-Meyer expansion fans which oscillate as 
the airfoil oscillates. These introduce serious complications into the 
second-order analysis. However, the complications can all be circum- 
vented by solving a "smoothed" problem in place of the actual problem. 
The solution can thereupon be applied to the actual problem, for which 
it yields the correct result everywhere except near the shock waves and 
Prandtl-Meyer fans. 



Sketch (b).- Steady flow past a 
wedge . 


The nature of the difficulties 
can be understood by considering 
first the special problem of steady 
flow past a single wedge (sketch (b)). 
The presence of the bow shock wave 
means that the analysis must be 
undesirably complicated by including 
the Rankine-Hugoniot relations (in a 
simplified form). A second compli- 
cation arises in the differential 
equation which, for steady flow, 
becomes 


t 


yy 


P^xx “ M 2 


(7+1) ft*®** + (7-1) ®x®yy + 24>y4>xy 


(15) 


In the iteration procedure to be employed, the nonlinear right-hand side 
is evaluated in terms of the first-order solution, and the resulting 
nonhomogeneous wave equation solved for the second-order potential. 
However, for the wedge the right-hand side vanishes (to any order), which 
would imply incorrectly that the second-order solution does not involve 
the adiabatic exponent 7. More precisely, the right-hand side vanishes 
everywhere except along the Mach lines springing from the apex, where it 
has the singular behavior of the Dirac delta function, and only by 
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taking account of these troublesome singularities could the correct 
solution be found. 

Both these complications are avoided by the simple device of 
solving the problem of flow past a smooth cusp-nosed airfoil of arbi- 
trary shape and then applying the final solution to the wedge. It may 
be imagined that the wedge has been 
smoothed by adding a cusped exten- 
sion to its nose, as indicated in 
sketch (c) . It is clear that this 
artifice removes the troublesome 
singularities from the right-hand 
side of the differential equation. 

Likewise, it eliminates the need 
for the shock-wave relations be- 
cause, as indicated in sketch (c), 
with sufficient smoothing, shock 
waves will form only at such great 
distance that their effects cannot 
reach the airfoil surface. Although Sketch (e).- Steady flow past 
shock waves are thus apparently a smoothed wedge, 

excluded, the correct second-order 

result for the wedge is nevertheless recovered from the solution by 
imagining the extension to shrink in size and disappear. The reason is 
that to second order a shock wave is equivalent to the limit of a rapid 
continuous isentropic compression. This limiting procedure, which is 
equivalent simply to applying the solution for an arbitrary smooth shape 
to one which is not smooth, yields the proper result except in the vicin- 
ity of the shock wave ( see ref. 5) • For an airfoil of general shape, 
similar broad smoothing must be imagined at any concave corner,* whereas 
at convex corners (since no shock waves form) the slightest rounding is 
enough. This smoothing technique was applied in reference 4 (p. 399) to 
steady first-order flow past bodies of revolution, and in reference 5 to 
steady second-order plane flow. 

We turn now to the question of generalizing this smoothing scheme 
to an oscillating airfoil. Modification is necessary only at the lead- 
ing edge. Consider first the special case of rotation about the leading 
edge. Then it is enough to conceive of an extension which is flexible, 
so that its cusped tip can be maintained fixed and directed always into 
the free stream while the airfoil 
oscillates, as indicated in 
sketch (d) . (The exact motion of 
the flexible tip is immaterial, 
provided the surface is suffi- 
ciently smooth and its slope 
remains small.) After the solution 
has been found, the flexible 
extension is again imagined to 




Sketch (d).- Smoothing for airfoil 
oscillating about leading edge. 
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shrink away, and the correct result is recovered for the act ual airfoil 
oscillating about its nose. 


y 



Sketch (e).- Smoothing for airfoil 
oscillating about arbitrary 
point . 


Finally, consider rotation 
about an arbitrary point. The flex- 
ible extension must now oscillate in 
such a way that its tip is always 
directed into the relative wind. 
Hence, an indicated in sketch (e), 
the tip must lie parallel with the 
free stream at the top (and bottom) 
of each stroke but incline in the 
direction of motion for intermediate 
positions. 


We are accordingly led to consider the motion of an arbitrary flex- 
ible oscillating surface described by 

y = e g(x) - ctoe iwt h(x) (l6a) 


where for the smoothed problem the functions g(x) and h(x) have con- 
tinuous first derivatives. The smoothed problem will ultimately be 
replaced by the actual problem. According to comparison with equa- 
tion (3b), this means that the function h(x) will eventually be identi- 
fied with (x - b). The requirement that the leading edge of the smoothed 
shape be always parallel to the relative wind may be written as 


g'(0) = 0 
h'(0) = iwb/U 


( 16 b) 


(TheBe last conditions, as well as the requirement that g' (x) be con- 
tinuous, must be relaxed in recovering the solution of the actual 
problem.) 


Boundary Conditions 


The boundary condition at the surface of the airfoil is that the 
normal component of velocity is zero. For any surface described by 
s(x,y,t) = 0 moving through a velocity field this condition means 
that the substantial derivative of S (i.e., its time rate of change 
for an observer moving with the fluid) vanishes at the surface (see 
ref. 6), so that 


~ 5 > 

St + q • grad S = 0 


at S = 0 


(17) 
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With velocity components given in terms of ® by equation (9b), and 
for the smoothed surface described by equation (l6a), this t ang ency 
condition becomes 

<I V = (! + ®x)( 6 S’ - c^e ^ h*) - ~ aoe h at y = eg - aoe iwt h (l8) 

where g s g(x) , etc. It is convenient to refer this condition to the 
axis y = 0 by expandin g in Taylor series . Keeping only terms of second" 
order gives 

®y = (l + 4 x )(eg’ 


(Here and 
ized theory.) 

The upstream boundary condition requires that in the actual problem, 
the Rankine-Hugoniot relations (or at least a simplified second-order 
form thereof) be satisfied across an oscillating bow shock wave whose 
position must be determined. However, shock waves have been eliminated 
from the smoothed problem, so that it is only necessary to require that 
the perturbation potential <S> vanish along the oscillating character- 
istic line (Mach line) springing from the leading edge. This insures 
that all disturbances produced by the airfoil are swept downstream. 

An equivalent and still simpler requirement is that $ and its stream- 
wise derivative vanish on, say, the plane x = 0: 

$ = <3> x = 0 at x = 0 (20) 


„ -iwt ico iwt , / iwt » 

- a^e h') - — c^e h - (eg - a^e h)®yy 

at y = 0 (19) 

on the right-hand side can be evaluated from linear- 


Restriction to Small Angles of Attack 


It is convenient to seek a solution to second order in the airfoil 
thickness, but to only first order in angle of attack. This is not a 
serious limitation because, just as in steady flow, second-order terms 
in angle of attack, though affecting local pressures, have no effect 
upon lift or moment because they are equal on the upper and lower 
surfaces. 

As a consequence of this restriction to linear terms in angle of 
attack, the awkward question of when the bow shock wave detaches will 
not arise. The amplitude of oscillation may be regarded as infinitesi- 
mal, so that if the bow shock is attached in the mean steady flow, it 
will remain so during the oscillation. 
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Transformation of Perturbation Equation 


It is convenient to separate the time -dependent part of the problem 
from the mean steady flow at zero angle of attack (for which the second- 
order solution is known) . Furthermore, for harmonic oscillations the 
number of independent variables is then reduced to two by separating an 
exponential time factor. Finally, the linear portion of the time- 
dependent equation is reduced to normal form by a transformation of 
dependent variable. These three transformations amount to setting 


®(x,y,t) = <t>(x,y) + Ooe 1 ^" 14 ^ fl) (x,y) 


where 


M 2 w 


(21a) 


(21b) 


Here <t> corresponds to the mean steady flow, and the term in $ repre- 
sents the additional flow associated with the oscillation through small 
angle of attack. 1 

Introducing this transformation into the perturbation equation 
(eq. (ll)) gives for 

*yy " P 2 *xx = M 2 ( 7 - + (♦/ + 4>y 2 ) x (22) 

where A is the Laplacian operation c^/ck 2 + ^/dy 2 , and for $ 

^ - (fj « * 

M 2 ^(7-l)( < t x A 1 ^+ $ X A4-) + 2 (<^X + 4^$y)x " iKl ^ (7-!)( 2 ^x + £►♦> + 

+ - k2 [2 + (r-1) M 2 ] * x 4 (23) 


•“•If second-order terms in angle of attack were retained, equation (21a) 
would include also a term c^ 2 e 2i(wt-Kx) X (x,y), which would appear 
only In the second-order solution. The differential equation for X 
would have the form 


X yy - - (Nr) 


x = }/F 


-1) t x A i|r+ 


and could be solved just as Is the equation for 
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The tangency condition of equation (19) likewise separates into the 
two conditions 


*7 = e (l + * X )g' - £ *yy S 


at y = 0 


(2k) 


= -e iKX (l + 4> x ) h 1 - iK iL e ittx h + e($ x - i^f) g* - 

e $yy g + e iKx 4>yy h at y = 0 (25) 

For the actual problem the second of these becomes, identifying h(x) 
with (x - b), 

= -e ilcX (l + 4> x ) - iK e lKx (x - b) + e (\j> -iK$) g* - 
J i/P x 


e$yy g + e iKX 4>yy (x - b) at y = 0 


(26) 


For pressures at the surface of the airfoil, the relation of equation ( ih-) 
can be expressed in terms of values at y = 0 by Taylor series expansion, 
with the result that to second order in thickness and first order in angle 
of attack 


C Pu 


(“2<t>x “ 2e<t> X y 
2o to e 1 ( UJt - KX ) 


g + p 2 ^ 2 - 




g + 


P 2<t, xfx ~ 4>y% + 




xy 



( 27 ) 


Solution by Iteration 


Although the equation for <t> is nonlinear, that for if is linear, 
but with nonconstant coefficients depending upon ♦. This corresponds 
to the physical concept that because of the restriction to linear terms 
in angle of attack, the oscillatory part of the flow is an acoustic 
field with, however, the speed of sound varying from point to point in 
accordance with the mean steady flow. 

The well-known linearized or first-order theory results from disregard- 
ing the right-hand sides of equations (22) and (23) . Thus, with the 
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first-order potentials denoted by the lower case letters cp and ty, the 
perturbation equations become 

<Pyy - ? ^XX = ° ( 28 ) 


%y “ P 2 ^xx 



(29) 


The second-order solution is obtained by iterating upon the first- 
order results. Using the linear equations to simplify the right-hand 
sides gives for the second-order iteration equations 


*yy - P 2 *xx = 2JF 5 [P 2 (N - l)<Px 2 + cp/] x (30) 

and 

^yy “ P 2 $xx “ ^ = ^ ^ (N - l) cpxtx + Ty'KyJx “ 

2iK [p 2 (2N-l)cp x ^ x+ P 2 NT xx f+cpyty]- 2N (jj)^ 


where 


B = 2L+J:g 

2 p 2 


(3D 

(32) 


(Here, following the usual subscript notation for derivatives, [ <P y 2 ] x 
means 8cpy 2 /8x, etc.) The second-order solution for 4>, which leads to 
Busemann’s well-known result at the airfoil surface (ref. l), was given 
in reference 5* It is therefore necessary to consider only the second- 
order problem for if. Details of the iteration procedure and discussion 
of its limitations are given for "the steady flow in reference 5 and 
apply also to the present problem. 


Partial Particular Integral 


The solution of the differential equation for steady flow in 
reference 5 'was simplified by discovery of a particular integral of the 
iteration equation in terms of the first-order solution. It was also 
shown there that for steady three-dimensional flow, a particular integral 
can be found to account for all terms in the iteration equation except 
those involving the adiabatic exponent 7 in the form of N. Likewise, 



NACA TN 2982 


13 


here, a partial particular integral which accounts for all terms on the 
right-hand side of equation (31) except those involving N is given by 

= & (<Pt) x - iKcpt (33) 

The complete- solution is this partial particular integral plus a 
solution of the reduced equation whose right-hand side contains only 
the terms still unaccounted for: 

*yy - P 2 *xx " 0[J $ = 2 P 2 n |i^(cfbctx) x - iK[(q>xt) x + <p x * x ] - <p x *} 

m 

First-Order Solution 


The first-order solution for <P is known from Ackeret's theory 

to be 

<P = - -j- g(x-Py) (35) 

It is to be understood here and in all similar expressions to follow 
that this is the potential only for x > py, and that cp vanishes 
identically ahead of the bow Mach wave (where x _<py). 

The first-order equation for ^ (eq. (29)) is most readily solved 
by applying the Laplace transformation with respect to x. We denote 
the Laplace transform of a function either by a bar, or by the symbol L, 
whichever is more convenient (and the inverse transform by L 1 ), so that 



Applying this transformation to equation (29), using the fact that 
and ijf x vanish at x = 0, gives 




( 37 ) 
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The solution of this equation which represents waves moving downstream 
is 


? = c(s) e nW s 2 + k2 /^ 


The constant C(s) is determined by the first-order form of the tangency 
condition (eq. (25)), which transforms to 


ty = - v(s) at y = 0 


(39a) 


where 


w(x) = e 1 * 1 * h'(x) + ik ^ h(x)J 


(39b) 


is the downwash velocity at y - 0. 
Consequently, 


_ 1 e-pyv^ 2 + K 2 /^ 

-»}r = — w(s) — 

P Jb* + K 2 /^ 


The inverse transformation is readily carried out using the standard 
tables (e.g., ref. 7) together with the convolution theorem, which gives 
as the solution of the smoothed problem 

x-Py 

* (x,y) = i J 3 o ^ V^-P 2 ;/ 2 J w(x -0 d| (1+1) 


The solution for the actual problem is now obtained by setting 
h(x) = x - b, which gives finally 2 

x-Py 

t(x,y) = -y J' J Q f i^ 2 -? 2 ^') e iK ( x- ^ ) 1 + Ik £ (x-|-b)J d£ (1+2) 

O \ s L- 

in agreement with the known result of linearized theory (see, e.g., ref. 8) 


The smoothing was, of course, unnecessary in this first-order problem. 
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SECOND -OEDER SOLUTION FOR LOW FREQUENCIES 


Because the first-order solution for arbitrary frequencies is 
rather complicated, use is sometimes found for an expansion in powers 
of frequency, which involves only elementary functions. The corre- 
sponding second -order solution will now be carried out in detail, 
including linear terms in frequency. This result will serve, for 
example, to evaluate the effects of thickness upon one-degree-of- 
freedom torsional instability, which is primarily a low-frequency 
effect. Thereafter, the result of extending the solution to include 
third powers of frequency will simply be stated. 


Potential Including Linear Terms in Frequency 


Expanding the first-order solution of equation (4l) in powerB of 
the frequency parameter k and retaining only linear terms gives 

z 

♦ (*,y) “if "(5) as + . . . o(k 2 ) = >Kz) + . . . (43) 

13 o 

where z = x - fty. To this order the partial particular integral of 
equation ( 33 )’ is a solution of the homogeneous equation (eq. ( 29 )), and 
can therefore be disregarded. Substituting the first-order solutions 
into the right-hand side of equation (3*0 and applying the Laplace 
transformation gives, to order K, 


^yy - 0 2 s 2 $ = -2eNe P 0 y (l^s-iK )L -jg' (x)ijr * (x)j^ -Ik s L -jg' (x)ijr(x)j- (*^0 


It is readily found that a particular integral of this equation is 
given by an appropriate multiple of ye“P s y. Then, adding a comple- 
mentary function representing downgoing waves gives 


$ = C(s)e-P s y + ye-P 8 y 


(f 



(*J-5) 


where the constant C is> to be evaluated from the tangency condition. 
Inverting the Laplace transformation shows that for the actual problem, 
in which h(x) = x - b, the solution has the form 


$ = f(z) + 



tf 2 g' + iK( 20 2 zg' - 0 2 bg r - g) 


(*t 6 a) 
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where g = g(z). The arbitrary function f(z) is determined from the 
tangency condition of equation (26) to he 


z , iK f P 2 ^ 

f(z) " s T (-3F* z -W h J 


+ e 


g + (*-„) 




z 

| 2-(W-l) / g(6)d| 


2(N-1) - 


21^-1 

iFp 5 J 


zg + z(z-b)g» + 


_1_ 

Iff 2 


(N-l) 



(46b) 


Surface Pressure Coefficient 


The pressure coefficient at the upper surface of the actual airfoil 
is found from equation (27) to be 


C Pu " C Po + 2o o e :L0t j- *) “ 


l^N-2 


eg' + 


Ik 


W 


_ N-l (2-l^)(#N-l) 

2M 2 —^2 g + £5 xg' + (l^N-2)bg* 


P 2 


♦ . . 0(K 2 ,e 2 a,a 2 ) 


(47) 


where g = g(x) . Here Cp Q is the value for the mean steady flow (at 
zero angle of attack) which is given by Busemann’s second-order theory. 
A more useful form of the result is obtained by extracting the real part 
and expressing the result in terms of the instantaneous angle of attack 
a(t) and its time rate of change a(t). Furthermore, the parameter e 
has served its purpose of distinguishing terms of different orders and 
can be eliminated (according to eq. (l)). Thus, on the upper surface 
of an arbitrary airfoil which is described at zero angle of attack by 
y = Y(x), is pivoted about a point a distance b downstream of its 
leading edge, and performs slow angular oscillations described by <x(t), 
the pressure coefficient is, to second order in thickness and first 
order in angle of attack. 
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°p u = c Po 



+ 


2 

0 



x + 



a 

U 


2 


lfN-2 

~W~ 


Y’a + 


2 


2^(N-1) Y + (2-MP)(lfN-l) xY , 
0 4 0 4 


+ !£2^bY' 
P 2 



Here the value for the mean steady flow is (ref. l) 



+ #^ Y ,a 

f ! 2 


0(a,e 2 a,a?) 

(48a) 


(48b) 


(a" preliminary report of this result was given in ref. 9„) In this 
form, the result is not restricted to sinusoidal motion hut applies to 
any oscillation which is sufficiently smooth and slow that the pressures 
depend significantly only upon the instantaneous angle of attack and 
angular velocity. 

The pressure on the lower surface of the airfoil is obtained from 
these equations by reversing the sign of a, and taking Y(x) to be 
the ordinate of the lower surface, measured positive downward. 

The result for plunging motions can be extracted by letting a 
tend to zero and b tend to infinity in such a way that their product 
remains finite, say 

ba(t) = h(t) (49) 

In the limit, the airfoil simply translates vertically according to 
y = h(t) . The pressure coefficient on the upper surface is 3 


C P 


u 


= Cp + a fi + 2 ififca 
•^o p u p 2 


Y 1 — 
U 


(50) 


Checks on the Result 


The solution can be teBted in several special cases for which the 
result can be derived from other considerations. 

3 Note that, as it should be, this is just the result of using Busemann's 
formula for steady flow (eq. (48b)), with the local slope increased by 

the instantaneous apparent downwash angle h/U; see the discussion in 
the following section of conditions at the nose. 
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Of the five terms in equation (48a), the first is known from 
Busemann’s steady second-order solution, and the second and third from 
linearized unsteady theory. The fourth is obtained by using the instan- 
taneous airfoil slope (T* -a) instead of the mean steady slope T* in 
Busemann's formula and retaining only linear terms in a. Therefore, 
only the last term, which is the essentially new result of the present 
analysis, requires verification. 

Just at the nose of an oscillating airfoil, the pressure can be 
determined exactly if the bow shock wave is attached. The transition 
through the moving bow shock is instantaneous, and so depends only upon 
the relative velocity at that instant (see ref. 4, p. 297). Hence the 
pressure just at the nose is instantaneously the same as on a wedge of 
the same vertex angle in steady flow with the same relative velocity. 

In the present problem, the relative 
velocity is compounded of the hori- 
zontal velocity U of the free 
stream, plus the instantaneous 
vertical velocity of the leading 
edge, which is given by db (see 
sketch (f)). The effect of the 
vertical component upon the equiva- 
lent free -stream velocity and Mach 
number is of second order in angle 
of attack, but the equivalent 
vertex angle of the airfoil is 
increased by the apparent downwash 
angle db/tJ. Replacing T* by Y* + abfu in Busemann's formula (eq. (48b)) 
gives, to first order in angle of attack 


U 



Sketch (f).- Velocity relative to 
leading edge of oscillating 
airfoil. 


] Pu ~ C Po 


2 ^ 
0 U 


+ 2 


l^N-2 


W 2 


(51) 


which checks the part proportional 
tion (48a) . 



Sketch (g).- Wedge oscillating 
about its vertex. 


to b of the last term in equa- 

The remainder of the term in 
question can be checked for a single- 
wedge airfoil oscillating about its 
vertex (sketch (g)). It can be shown 
using the results of reference 3 
that in this case disturbances 
reflected from the shock wave are of 
third order in the wedge angle 
(although for other pivot positions 
they are of second order) . There- 
fore, a solution correct to second 
order in thickness and first order 
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in angle of attack can "be found by applying linearized theory to the 
mean steady flow behind the shock wave. 4 For slow oscillations, the 
first three terms of equation (48a) give 


P-Pi 


2 


= a + 

01 


2 (2 -Ml 2 ) 

0l 3 


Xl 
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Ui 


(52) 


where subscript 1 denotes values in the mean steady flow behind the 
shock wave. From linearized theory 
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(53) 


where e is the semi vertex angle. Hence, referring the surface pressure 
coefficient to free -stream quantities (and noting that x x - x to second 
order) gives 



- §• a +2 
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2-M 2 

0 3 



- 2 


KpN-2 

ea + 

0 2 


„ 21^(W-l) + (2-l^)(^N-l) 
2 13 



( 5*0 


which checks the first two parts of the last term in equation (48a) 
when Y(x) = ex. 

Recently, Lighthill has given a further check for the case of 
Mach numbers so high that l/tf is negligible compared with unity 
(ref. 10) . 

This concept was suggested to the author by W. P. Jones of the National 
Physical laboratory, England. 
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Lift anri Moment Coefficients for Symmetrical Airfoils 


The coefficients of lift and pitching moment (about 
given in terms of the pressure coefficients on the upper 
faces by 

c 



the pivot) are 
and lower sur- 


(55a) 


1 





(55b) 


For simplicity, consider only airfoils symmetric about the chord line. 
In this case, the pressure difference is given by 


C PZ - C Pu 


k- 1 V f 2-M 2 h \ d 


+ Y’a- 

P 2 


2*f(N-l) Y + (2-^)(tf»K-l) ^ , + M*N-2 bY , 

- p 4 P 4 P 2 


- (56) 
U 


If the airfoil has a sharp trailing edge, substituting into equations 
(55) and integrating by parts gives 
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(57) 
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/ Ydx 
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2-M 5 , 2M 2 -3 b 
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b _ b£\ _ 

c c z J 


ca 


(58) 


Inside each square bracket the first term is the result of linearized 
theory, and the remainder represents the second-order effect of thick- 
ness. Thickness effects are seen to appear in tne form of the area of 
the airfoil profile and its first moment about the vertical line through 
the pivot. 
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If the airfoil has a blunt trailing edge of semi thickness Y(c), 
the following additional, second-order terms must be added to the above 
expressions: 


to cj: 


k M 2 N-2 Y(c) _ k 

P P c a (3 


t^N-2 b . (2-M 2 )(l^N-l) 
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(59) 


to cm: 

k f^N -2 b\ y(cj_ a + 4 
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M^N-l) 

P 3 


I^N-2 


_ b \ “1 4 b\ Y(c) cd 

c J J \ cj c U 


(60) 


Examples: Biconvex and Double-Wedge Airfoils 


To second order, a biconvex airfoil of thickness ratio t is 
given by the parabolas 

y = ±Y(x), Y(x) = 2 t ~ ( c -x) (6l) 

The expressions for lift and moment become 
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(63) 


It happens that for a double-wedge airfoil (with maximum thickness at 
midchord), both the area and first moment are just three fourths of 
those for the biconvex airfoil. Consequently, the above results apply 
to double-wedge airfoils, if t/3 is replaced throughout by r/b. 


In the expression for pitching moment, the term proportional to a 
represents an aerodynamic stiffness or restoring moment in phase with 
the angular displacement, while the term proportional to d corresponds 
to an aerodynamic damping moment in phase with the angular velocity. 

The effects of thickness upon aerodynamic restoring and damping moments 
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are shown in sketches (h) and (i) for a 6-percent-thick double-wedge or 
a 4-l/2-percent-thick biconvex airfoil with two different pivot posi- 
tions. The sketches have been labelled with both the usual American 
notation (e.g., ref. 11) and the British notation (e.g., ref. 12), 



Li nearbed theory 

Second- order theory 


1 

. 

“ 


b/c * 

’^oli ST 

s 

A B / L / 

Zi N 6 


i.OO LE5 iso L75 2.00 

M 


Sketch (h).- Effect of thickness upon restoring moment coefficient for 
4- l/2 -percent-thick biconvex or 6-percent-thick double-wedge 
airfoil . 
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Linearized theory 

Second-order theory 
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* e.5 
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A B / L / 
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Sketch (i).- Effect of thickness upon damping moment coefficient for 
I 4 .-I/ 2 -percent-thick biconvex or 6-percent thick double-wedge air- 
foil. 


defined for harmonic oscillations of arbitrary frequency by 
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Here k and A are the reduced frequency in the American and British 
notations, respectively, related hy 


A 


= 2k = 


OJC 

U 


(65) 


For slow oscillations, 


the coefficients are given "by 
2 oa 


k - _ U 8c m 
md = _ Mi = ^ 


Sd 


(66) 


(and, indeed, this was originally the definition of and m<x, that 
given by equation (64) being a later extension to the case of rapid 
harmonic oscillations). 


It should be noted that according to second-order theory, the non- 
linear effects of thickness are themselves linear in thickness. This 
means, for example, that doubling the airfoil thickness ratio would 
double the distance between the linearized and second-order curves of 
sketches (h) and (i). 


Neutral Damping Boundary 


Linearized theory indicates the possibility of instability of 
pitching oscillations for low frequencies. For a range of Mach numbers 

below a/5/2 = I .58 and pivots ahead of two thirds of the chord, the 
aerodynamic damping moment becomes negative, and so tends to destabilize. 
(Whether or not the motion is actually unstable depends, of course, upon 
the other dynamic parameters in the problem.) This zone of possible 
instability shrinks and eventually disappears as the frequency of oscil- 
lation increases. The present low-frequency solution is therefore 
adequate for determining how the region of instability Is modified by 
nonlinear thickness effects. 

Sketch (j) shows the boundary of neutral aerodynamic damping for 
slow oscillations of' a 4- l/2-percent-thick biconvex or 6-percent-thick 
double-wedge airfoil. The aerodynamic damping is destabilizing for 
Mach numbers and pivot positions lying inside the loops. Within the 
region where linearized theory predicts a destabilizing moment, thick- 
ness is seen to exert a further destabilizing effect except for pivots 
near midchord. The second -order solution becomes unreliable when the 
bow shock wave detaches, at about M = 1.2. 
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Sketch ( j) . - Neutral damping boundary for 4-1/2 -percent -thick biconvex 
or 6-percent-thick double-wedge airfoil. 


Comparison With Previous Investigations 


Two previous investigators have sought a second-order solution for 
slowly oscillating airfoils in supersonic flow. Their results agree 
neither with each other nor with the present solution. 

In 1947, W. P. Jones obtained an estimate of the thickness effect 
by assuming that the ratio of second-order to linearized pressure 
disturbances is the same for slow oscillations as that given by 
Busemann's formula for steady flow (ref. 13) • That this assumption is 
not altogether correct is indicated by the fact that the results do not 
check those obtained for a wedge oscillating about its vertex by apply- 
ing linearized theory to the mean steady flow behind the shock wave. 
However, the assumption is correct at the leading edge, and also (as 
noted by Lighthill in ref. 10) in the limit of high Mach number. It is 
seen in sketch («j) that thiB estimate fails to give a useful prediction 
of the actual effects of thickness, except for pivots near midchord. 

In 1951 j Alexander Wylly attacked the problem by methods similar 
to those used here (ref. l4) . Unfortunately, it appears that the 
smoothing was not carried out with sufficient care; as a consequence. 
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the solution satisfies none of the three checks discussed previously. 
In contrast to the present results, the effect of thickness upon aero- 
dynamic damping was predicted to he stabilizing an d so great that for 
airfoils of the thicknesses shown in sketch (j), the zone of possible 
instability would have disappeared altogether. 


Comparison With Experiment 


The only experimental results available are those quoted in 
reference 15, which are derived from tests by Brairt and Chinneck. 

A 7-l/2-percent-thick biconvex airfoil, pivoted at midchord, was tested 
at Mach numbers of 1.28, 1.46, and 1.52 over a range of frequency. The 
aerodynamic damping was measured by giving the airfoil an initial dis- 
placement and measuring the logarithmic decrement. The experimental 
data are shown in sketch (k) in comparison with linearized theory and 
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L2 


1.3 

¥ 


1.4 


1.5 


L6 


Sketch (k).- Comparison of experimental and 
theoretical damping moment coefficients for 
7-l/2-percent-thick biconvex airfoil pivoted 
at midchord. 

the second-order theories of W. P. Jones, Wylly, and the present 
analysis. As indicated on the sketch, the lowest test Mach number lies 
in the range of mixed subsonic and supersonic flow (and a detached bow 
wave). The present theory predicts only a fraction of the discrepancy 
between experiment and linearized theory. The same is true of Jones’ 
estimate which, according to sketch (j), agrees with the present analysis 
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nn t y for pivots near midchord. Wylly’s result, which has been shown to 
be erroneous, happens to lie much closer to the measured results. 


Extension to Cube of Frequency 


The dependence of thickness effects upon frequency of oscillation 
can be est ima ted by extending the second-order solution to include 
higher powers of frequency. This has been carried out for an arbitrary 
airfoil by incl uding second and third powers, which is enough to show 
an effect of frequency upon both aerodynamic stiffness and damping. 

The computation, though cumbersome, is a straightforward extension of 
the previous analysis, so that only the final result will be given here. 


The expression for pressure coefficient on the upper surface, 
corresponding to equation (4-7) , is found to be 
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where g = g(x) and / g = /g(£)dg, etc. (The first four terms are 

o o 

the result of linearized theory.) This result meets the three tests 
discussed previously, arid also checks the solution given later for a 
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vedge with general pivot location at arbitrary frequency. The result- 
ing expressions for lift and moment involve the airfoil thickness in 
the form of the area of the profile and its first three moments about 
the vertical line through the pivot. 


Example: Biconvex Airfoil 


These rather formidable results simplify considerably for specific 
airfoils. For example, for a biconvex airfoil of thickness ratio r 
oscillating about midchord, the pitching-moment coefficient is given by 


c m 
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The first three terms are the result of linearized theory. 

The component of this moment which is out of phase with the angle 
of attack gives a parabolic approximation for the variation of aero- 
dynamic damping with frequency as shown by dashed lines in sketch (Z) 



O .1 .2 .3 .4 .5 .6 .7 

Reduced frequency. A* 2k * *c/U 

Sketch (Z).- Effect of frequency upon damping moment coefficient 
for 5 -percent -thick biconvex airfoil pivoted at midchord. 
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for a 5 -percent -thick biconvex airfoil. The accuracy of the parabolic 
approximation for linearized theory is indicated by comparing it with 
the exact result (solid line). In this example the linearized and 
second-order curves run almost parallel, which means that the nonlinear 
effects of thickness vary only slightly with frequency. 

Recently, Jones and Skan have treated biconvex airfoils at arbi- 
trary frequency by a numerical procedure (ref. 16) . Their result is 
shown in sketch ( l ) . It fails to give the initially parabolic form of 
the curve which is implied by the fact that the second-order solution, 
like the first-order result, can be expanded in powers of the square of 
the frequency. Their solution involves several doubtful assumptions, in 
particular, that the effect of the bow shock wave can be disregarded. 

It has already been remarked here that, actually, the bow shock has a 
second-order effect unless the pivot lies at the leading edge. 


SECOND-ORDER SOLUTION FOR ARBITRARY FREQUENCIES 


For some purposes the previous solution for slow oscillations may 
be inadequate. In principle, the second-order solution can be extended 
to include still higher powers of frequency, but the labor required is 
clearly prohibitive. Alternatively, one can attack directly the problem 
for arbitrary frequencies. 

The second-order solution can, in fact, be carried out for a general 
airfoil at arbitrary frequencies. However, the result is formidable, 
involving multiple integrals of products of Bessel functions, and the 
reduction to simpler form for specific profiles appears to be difficult. 

A more practical approach is to choose a specific airfoil shape in 
advance. Then the second-order solution involves only functions of the 
type encountered in the linearized theory. In particular, it is found 
that (at least for the simplest shape) the final expressions for lift 
apfl moment involve only functions which have been already studied and 
tabulated, so that numerical results are readily obtained. 


Modified Smoothing Procedure 


The smoothing discussed previously must be dropped at an earlier 
stage of the solution when a specific airfoil shape is chosen. It is 
therefore necessary first to modify the differential equation and 
boundary conditions so that no singular terms appear. 

Consider first the differential equation (eq. (31)). Applying the 
Laplace transformation of equation (38), and envisioning the smoothed 
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problem, so that “$■ and vanish at x = 0, reduces it to 

fyy ~f (S 2 + ^ f = 2(^S-i K ) (N-DX^^j. + I^Tyj - iK p j] 

(69) 


Tn this form, an troublesome second derivatives have disappeared from 
the right-hand side, so that the smoothing can now be dropped. 

Consider next the tangency condition of equation (25) , which 
contains second derivatives of both cp and These cause no difficulty 
when, as in the solution for slow oscillations, the tangency condition 
i's imposed as it stands. However, its Laplace transform will be used 
here, and then the Dirac delta functions associated with the second 
derivatives would affect the integration implied in the inversion of the 
transformation. These troublesome second derivatives can be eliminated 
by first expressing them in terms of x derivatives only through the 
first-order equations: 


^yy " P 2 ^cx 

V * <♦= + § ♦) 


(70) 


and then applying the Laplace transformation while envisioning the 
smoothed problem (so that <p, and i|r x vanish at x = 0), which 

gives 

- w(s) - + e^Irjfeg'j- - lice L-jtg’j- - pPes^tacsJ- - 

P 2 ^ L-jtgl + P 2 ( s-iK ) L-fe iKx cp x (x-b)j- at y = 0 (71) 


Again, all second derivatives have disappeared, so that the smoothing 
can-be dropped. 

Consider finally the upBtream condition. For the smoothed problem, 
one statement of the condition was seen to be that the solution repre- 
sents dawngoing waves. ThiB means that the complementary function for 
the iteration equation should have the same form as the linearized 
solution, and this statement of the upstream condition applies as well 
to the actual problem. 
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Solution for Wedge 


The simple case of a single wedge illustrates the method of solu- 
tion which can, in principle, he applied to any profile formed of 
piecewise analytic arcs. For a wedge of semi vert ex angle e, the 
first-order solution for cp is given by equation (35) as 

q> = - 5 (x - Py) (72) 

3 

and the first-order solution for '4 r is given by equation (42). Substi- 
tuting into the right-hand side of equation ( 69) 5 gives for the trans- 
formed iteration equation 
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A particular integral of this equation is given by an appropriate 

multiple of ye Py^/ 8 +K 2 /M 2 ^ and adding a complementary function which 
represents downgoing waves gives the solution 
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The constant C is evaluated by imposing the tangency condition of 
equation (71)# with the result that the Laplace transform of the second- 
order solution is found to be 
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^Cn this simple case, the partial particular integral of equation (33) 
offers little actual simplification « 
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The inversion can "be carried out using the standard tables (e.g., 
ref. 7) together -with the convolution theorem. For calculating surface 
pressures, it suffices to obtain the solution in the plane y = 0, 
which is found to be 


f (x,0) B J J o ^ ^ v(x-£)<U + ex v(x) + 

o 

e £ ^W°) J ° (I *) - J ° (5 *) + 
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j? f J o 0{ 1 V*(^)-tt(2N-l)w(x-S) - 


x-i 

K * i f w( 11)471 ] - eiK(x " 6) } 


(76a) 


where 


w(x) = e ilcx 


1 + iK 


h (x_b) 


(76b) 


With the pivot at the nose (b =0), this agrees with the result of apply- 
ing linearized theory to the mean steady flow behind the shock wave. 

Also, when expanded in powers of frequency, it agrees with the previous 
low-frequency solution up to terms in It 3 . 

The surface pressure coefficient can now be calculated from equa- 
tion (27) and the lift and moment coefficients from equations (55) • 


Example: Wedge Pivoted at Nose 


For simplicity, the results will be given only for the special 
case of rotation about the nose. Then it is found that the lift and 
moment coefficients are given by 
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where v = 0 for the lift and v = 1 for the moment. The functions 
Tel, given by 


fn(M, A) 


J x^-i^x/P 2 Jo 



(78) 


arise in linearized theory for n ranging from 0 to 3. They have been 
studied and tabulated by von Borbe'ly (ref. 17), Schwarz (ref. 18), and 
Garrick and Rubinow (ref. ll) . They can all be expressed in terms of 
f Q "by a recurrence relation due to von Borb^ly (see* e( 3L* (A*87) 

of ref. 12), so that the additional f 4 required here is easily 
computed. 
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Sketch (m) shows the variation of aerodynamic damping moment. 



O A .8 L2 L6 2.0 2.4 2.8 


Reduced frequency, X =2k* mc/u 


Sketch (m).- Damping moment coefficient for 10-percent- 
thick wedge oscillating about its vertex at 
M = 10/7 o 


according to first- and second-order theories, for a 10-percent-thick 
single wedge oscillating about its nose at a Mach number of 10/7. Also 
shown for comparison are the parabolic approximations of the low-frequency 
a n alysis. It is seen that the thickness effect is reversed at high 
flutter frequencies, as is suggested by the parabolic approximation. 


DISCUSSION 
Higher-Order Effects 


The moderate magnitude of second-order effects would suggest that 
the influence of third- and higher-order terms is of no practical 
importance, except perhaps in the transonic range near shock detachment. 
This supposition can be confirmed in the case of the single-wedge airfoil, 
for which a solution exact in thickness (but linearized with respect to 
angle of attack) has been derived in references 2 and 3* Sketch (n) 
compares the boundaries of neutral aerodynamic damping for a slowly 
osci llat ing wedge of 5° semivertex angle as predicted by the linearized. 
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Pivot position, b/c 


Sketch (n) Neutral damping boundary for wedge of 5° 

semi vertex angle. 


second-order, and exact theories. The second -order solution lies close 
to the exact result down to the Mach number for shock detachment (which 
is almost the same as the Mach number at which the flow ceaseB to be 
purely supersonic ) „ 


Application to Finite-Span Wings 


Extension of the second-order solution to wings of finite aspect 
ratio does not seem possible at present. No second-order solution has 
yet been found even for steady flow past the simplest lifting wing. 

Fortunately, the main conclusion to be drawn from the present 
analysis is that nonlinear thickness effects are quite moderate in 
magnitude. Practical supersonic wings will, therefore, probably be 
so thin that nonlinear effects are negligible, so that reliance can be 
placed in the predictions of linearized theory. Only if the wing is 
unduly thick, or if the Mach number is close to unity, or if unusual 
accuracy is required, may the engineer be forced to estimate the effects 
of thi ckness - In this event, he might assume that the effects of thick- 
ness are in some sense additive to those of aspect ratio, provided the 
aspect ratio is high and the frequency low. For example, the two- 
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dimensional correction might he applied stripwise to the spanwise load- 
ing predicted by linearized theory. Some indication of the extent to 
which such an assumption would be justified can be obtained by consider- 
ing other pairs of effects whose combined influence is known. Sketches 
(z) and (m) show that the effects of thickness and frequency are roughly 
additive for the frequencies of usual practical interest in dynamic 
stability analysis (say A<0.2). Likewise, sketch (o) shows that the 



Reduced frequency, X = 2k = *c/u 


Sketch (o).- Effect of aspect ratio upon damping moment, 
coefficient for rectangular wing according to 
linearized theory. 


effects of aspect ratio and frequency, determined from Watkins • linear- 
ized solution for the rectangular wing (ref. 19), are nearly additive 
in the same range of frequencies. 

Recently, Martin and Gerber have calculated the second-order effects 
of thickness upon damping in roll for wings of infinite span (ref. 20). 
They then estimate the effect for a finite rectangular wing by increasing 
the result of linearized theory in the same ratio as for the infinite 
wing. Agreement with experiment is thus considerably improved, which 
gives further assurance that superposition of thickness effects may be 
justified also for oscillating wings. 
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Further Analysis 


If thickness effects of the magnitude indicated hy the present 
analysis are judged to he significant at flutter frequencies, extension 
of the high-frequency solution to more practical profiles would he 
warranted. The solution given for a single wedge at arbitrary frequency 
should he extended next to the biconvex or double -wedge airfoil. 

Although considerable computation is involved, it does not appear that 
the labor would be prohibitive. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., April 20, 1953 
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APPENDIX 

SYMBOLS 


A 

b 

c 

C(s) 

c l 

c m 

°P 

C Po 

C Pu' C Pl 
g(x) 

b(t) 

h(x) 

k 

®a> Mj 

M 

N 

P 

q. 


speed of sound 
aspect ratio of wing 

downstream distance from leading edge to pivot 

airfoil chord 

constant of integration 

section lift coefficient 

section moment coefficient 

pressure coefficient 

surface pressure coefficient in mean steady flow 

pressure coefficients on upper and lower surfaces of 
airfoil, respectively 

function defining upper surface of airfoil at zero angle of 
attack 

elevation of plunging airfoil 

function defining amplitude of oscillation of upper surface 
of smooth airfoil 

toe 

reduced frequency, — 

2U 

aerodynamic damping coefficients 
(See eq.. (66)) 

free-stream Mach number 
7+1 

2 P 2 

static pressure 
velocity vector 

Laplace transformation variable 
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S(x,y,t) function defining moving surface 
t time 


u,v 

U 

w(x) 

x,y 

Y(x) 

z 

a 

«o 

P 

7 

A 

e 

K 


velocity components parallel and perpendicular to free stream 
free-stream velocity 
dovnwash velocity at y = 0 

coordinates parallel and perpendicular to free stream 
ordinate of upper surface of airfoil at zero angle of attack 
x - py 

angle of attack 

angular amplitude of oscillation 

y m 2 - 1 


adiabatic exponent of gas 

Laplacian operator, + 

Sy 2 

small parameter representative of airfoil thickness 


& 

p 2 U 


Tv reduced frequency, uc 

U 

v constant which is 0 or 1 in equation (7*0 

p density 

t airfoil thickness ratio 

9 first-order mean steady perturbation potential 

<t> second-order mean steady perturbation potential 

G> complete perturbation potential 

first-order time -dependent perturbation potential 
\J/ second-order time -dependent perturbation potential 
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partial particular integral of time-dependent iteration equation 
angnl ar frequency of oscillation 
complete velocity potential 
free- stream conditions 

conditions in mean steady flow “beh in d bow shock wave on wedge 
differentiation with respect to time 

Laplace transform 
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